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Abstract 

Trigonometric invariants are defined for eacli Weyl group orbit on the root lat- 
tice. They are real and periodic on the coroot lattice. Their polynomial algebra is 
spanned by a basis which is calculated by means of an algorithm. The invariants 
of the basis can be used as coordinates in any cell of the coroot space and lead to 
an exactly solvable model of Sutherland type. We apply this construction to the F4 
case. 
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1 Introduction 



Integrable models of the Calogero-Moser class and their trigonometric and rational limit 
models are conventionally described by the simple Lie algebras, i.e. those contained in 
the classical sequences y4„, C„, D„ or the exceptional set G2, F4, Eq, E7, Eg. We have 
shown 3 that the Weyl groups underlying these algebras are the essential ingredients 
in the construction of the trigonometric or rational models. Exact solvability is easily 
proven by expressing the Schrodinger equation in terms of the Weyl group invariants as 
coordinates. For the rational models the crystallographic property of the Weyl groups, 
which guarantees the existence of a root lattice, can be abandoned. For the proper Coxeter 
groups H3 and and the infinite sequence of dihedral groups 0, 0] Calogero type models 
can be constructed as well. 

For the classical Lie algebras and G2 the Weyl group invariants are not derived but 
guessed by intuition But this method failed for the exceptional algebras F4 and 
Eq, Ef, Eg. Therefore a concept supplying us with all trigonometric polynomial invariants 
and their algebraic basis is highly desirable. In the case of the rational polynomial invari- 
ants the existence of the algebraic basis is guaranteed by Chevalley's theorem 0. The 
Jacobian for the transition from cartesian coordinates in root space to the basic invari- 
ants as new coordinates can be factorized ("factorization theorem"). Both the Chevalley 
theorem and the factorization theorem are valid also in the trigonometric case as we shall 
show here for F4 (for the classical Lie algebras and G2 they are known to be valid, too). 

Our construction of invariants proceeds as follows. We decompose the root lattice 
A into orbits Q (infinitely many, they can be ordered by the length of vectors which is 
constant over the orbit). For each orbit we define an invariant trigonometric polynomial 




(LI) 



pen 



These functions obey fusion rules 




(1-2) 



with "fusion coefficients" C^^^ that are nonnegative integers. In ( |L2| ) the null-orbit 



consisting only of the null-vector in A must be included 



Tno(x) = 1 



(1-3) 



The system of equations (|1.2| ) is of triangular shape and can be solved trivially for 
^n,nax(^)- Each pair Qi,Q2 defines in fact a unique i7jnax(^i7 ^2) with 



(f2max(Sll,Sl2) 



1 



(1-4) 



By recursive substitutions we isolate then an algebraic basis 




(1-5) 



and for all other orbits we obtain explicitly 



Tn(x) = pol{Tni(a;),...TQ„(a;),TQo(x)} 



(L6) 



This is an explicit and constructive version of Chevalley's theorem. It is obtained only in 
a case-by-case study (e.g. for F4). In Section 2 we do this in great detail for F4. 

In Section 3 we apply this technique to construct the F4 Sutherland model, using the 
approach developed in 0]. By the way some minor theorems are proven by explicit 
calculation (e.g. existence of the rj-"^-coefficients as polynomials in the Chevalley basis). 

2 Trigonometric invariants of Weyl groups. 

Weyl groups are generated by reflections along roots a G ]R„ p 

X e Rn ■■ SaX = X - (2.1) 

(a, a) 

Each root a is an integer linear combination of simple roots 

{ai,a2, ...an} (2.2) 
These simple roots span an integral lattice A G 

n 

13 e A: p = ^miai, rrii e Z (2.3) 

i=l 

If the Weyl group W acts on a vector /3 G A it produces an orbit Q 

n = {w(3, weW} (2.4) 

How can such orbit be characterized? 

Since is a discrete subgroup of 0(n) acting on ]R„ we obtain 

= {wf3,w(3) 

= m' (2.5) 

so that an orbit appears as a discrete set on a sphere of radius and (see ( p.3|) ) 

II/^IP = X]*^"" (2-6) 

Here (ctj, aj) is contained in the Cartan matrix as 

A, = 2^^^ (2.7) 

Only for simply laced Weyl groups the lengths of the simple roots are all equal. 
For the non-simply laced Weyl group W = F4 we use as basis in R4 

{ci, i G {1,2,3,4}, {ei,ej)=Sij} 

Then W = can be generated from the reflections along Cj and fi (all i) where 



fl = 


1/ 

-(ei + 62 + 63 + 64) 




/2 = 


1/ ^ 
-(ei + 62-63- 64) 




/3 = 


-(ei -62 + 63- 64) 




h = 


-(ei -62-63 + 64) 


(2.8) 



if^Jj)=5^J (2.9) 

and the permutation group 5*4 of the basis {cj}. The set of roots decomposes into two 
orbits 

= {±e,, ^(±ei ± 62 ± 63 ± 64)} (2.10) 

f^2 = {±e.±e,, i<j} (2.11) 

which are characterized by a subscript denoting the length squared of the roots. In 
addition we need the null-orbit 

Qq = {null- vector} (2.12) 

All other orbits can be characterized by an integral radius squared and a further "degen- 
eracy label. A list of orbits up to ||/3|p = 24 is given in Table 1. 
Trigonometric invariants are defined for each orbit by 

Tn,a{x) = ^ expi{P,x) (2.13) 

(i.e. Wpf = n), so that 

T„,a(0) = (2.14) 
Now we consider pairs of such trigonometric invariants and expand their product as 

k,c 

Applying permutations of 5*4, reflections along the coordinate axis Cj, i G {1,2,3,4} and 
reflection along /i to any generating vector yields the whole orbit. 



We expand x G M„ in a co-root basis 



Then by inversion of (|2.15|) we obtain 



(27r)-" / rf"er„,a(x)r„,,(x)T,,,(x) 

Jcell 



Note that the r.h.s. is symmetric in the three orbits. We denote ( p.l5| ), (|2.17|) the "fusion 



rules" for trigonometric invariants. The "fusion coefficients" C|^^^^ are nonnegative 
integers. 

We expect that a Chevalley theorem of the following type is valid: 
Any polynomial in the {T„ a(a;)} can be expressed as a polynomial in an algebraic basis 
of invariants including Tq^ = 1. The number of nontrivial basis elements is rank W. 



In the case of W = F4 this algebraic basis is constructed by inversion of (|2.15|) and 
consists of 

T,,T,,Ts,TP (2.18) 

I This inversion is possible by the triangular shape of the fusion rules. Namely inserting 
( |2.13|) into ( |2.15| ) and using the triangular inequality we obtain 

C[nS){m,b) = except possibly for 

\\fn — \fm\ < < \fn ^ \fm (2.19) 

Thus there is a maximal k for each n, m 

(2.20) 

and we solve ( |2.15D for Ti^^^_^^J^x) where c is such that the fusion coefficient is one. The 
result for 

Tfc,„ fc^ {1,2,3,6}, A;<24 (2.21) 

expressed as a polynomial in Ti, T2, T3, Tg is given in Table 2. 

The task to introduce the trigonometric invariants as coordinates in each cell of the 
space R„ leads to the study of the Jacobian matrix 

f^l (2.22) 

* ) mgbasic set of invariants 

In the program of constructing exactly solvable models we have to compute the (inverse) 
Riemannian 

-1 _ \ ^ <9Tm, dTn , ^ 

n 

k,a 

which can obviously be expressed as a polynomial in the basic invariants. For F4 the 
largest orbit appearing is ^^24 in g^Q- 



*) We could replace Te by one of the following invariants: {T7, Tg, Tg^i, T9_2, T'lo} or constant linear 
combinations thereof. 



3 The F4 Sutherland model 



From ( [^.23D a first version of tlie Riemannian is obtained by insertion of tlie fusion coef- 
ficients 

911 = -T4 - T3 + 4Ti + 24 

912 = -^5 + 6T1 

9n = -^T7-3T6-2T5 + 3T3 + 12T2 + 12Ti 

-1 



(^re = -2Tn - 2T9,i + 4T5 + eTg 
g^i = -2T8 - 2T6 + 8T2 + 48 
g^^ = -2T9,i - 3T7 + 6T3 + 24Ti 
g^e' = -3Ti4 - 6T12 - 4rio + 6T6 + 24r4 + 24T2 
^3-3^ = -3Ti2 - 5Tn - 8T10 - 12T9,2 - 6X9,1 - 12T8 - 3T7 + 4T5 

+ I2T4 + 2IT3 + 48T2 + 6OT1 + 288 
^3-6^ = -4Ti7,i - eTis - 6Ti3,2 - 8ri3,i - 4Tn 

+6T7 + I6T5 + 36T3 + 48Ti 
9m = -6T24 - IOT22 - I6T20 - 12Ti8,2 - 24Ti8,i 

-24Ti6 - 6T14 + 8rio + 24T8 + 42r6 + 96T4 + I2OT2 + 576 (3. 1) 

A second form is obtained by substitution of the algebraic basis Ti, T2, T3, Tg by means of 
Table 2. 

= -Ti' + r3 + 6T2 + 12Ti + 48 



g-^^ = _TiT2 + 3T3 + 12Ti 

^1-3' = -^TiT3 + 4TiT2 + 12Ti2 + ^r6-24r3-42T2-96ri-288 

g^i = -2T1T6 + 2T2T3 + 4T1T2 - 12r3 - 48Ti 

c/22^ = -2T^ + l2Tl + 2T6 - 24T3 - 48T2 - 96Ti - 192 

^2-3' = -2T2r3 + 3T1T3 - 4T1T2 - 24r2 - 9T6 + 6OT3 



+ IO8T2 + 240Ti + 576 



926 



- -ATIT2 - 2ATl - 3T2TQ + 3T| + 8r3T2 - 6r6Ti 
+72T3T1 + IO4T2T1 + 12T| + 96T^ - 48T6 + 288T3 
+48OT2 + 1536Ti + 2304 
g^i = AT2T^ + 12T^-AT2T,-3T^ + nT^-3QnT, 

-6OT1T2 - 96T^ + 12Te - 48T3 - 48T2 - 384Ti 
^3-6^ = 2T1T2T3 - 16T2T2 - 4T3T6 - 6T2T6 - 8T1T6 + 72T| 
+32T2r3 - I2T1T3 + 144rir2 + 96T^ + 36r6 
-24OT3 - 48T2 - 960Ti - 2304 
g^i = -16T^T2 - mT2n - 8T^n + 2T2T^ 

+48T1T2T3 - 96T2T2 + 48TiT| - 192T^^ 
-6T2 + I6T3T6 - 24r2r6 + IQTiU + 2ATi 
+512T2T3 + 576T1T3 + 864T| + 288OT1T2 + 1344Ti^ 
+96T6 + II52T3 + 6I44T2 + 7680Ti + 9216 



(3.2) 



Since has two orbits in the roots we have according to the factorization theorem 



where Pj corresponds to Qi ( p.lO| , |2.11|) . We find exphcitely 



Pi = 110592Ti + 41472T2 + 27648T3 + 110592 - 3456T6 - 1728TiT6 
+192TiT| - 432T3T6 - 112Ti^T3 - SSAT^T^ 
+5I84T2T3 + 20736T1T2 - n2ST^T2 + UAT^Tq 
-48Tj^T2 - 648T2T6 + I4976T1T3 + ATqTI 
-TItI + 2I6T1T2T3 - ISTqTiTs - 4608Tf 
+1728T| + 18432T2 + 3888T| + 27T| + 16Tf 
+4T| 

P2 = 10616832Ti + 4423680r2 + 1769472T3 + 7077888 - 221184r6 
-221184TiT6 + 103680TiT| - 27648T3T6 
+34560T3T3 + 663552r2T3 + 774144T2T3 
+4866048Tir2 + 1465344X^2X2 - 62208X1^X6 
-6912XfX2 - 78336X2X5 + I99O656X1X3 
-1728X6Xf + 18144X|X2 - X^Xg^ 
+566784X1X2X3 - I728OX6X1X3 - 8X6XiX| 
+79488X1^X3X2 - 48384X6XiX2 - 2592X5X1X1 
+43200X1X1X3 + 36X5^X1X2 - 2592X3X1^X6 
-8640X3X1^X2 - 46O8X2X3X6 - 5184X2X1^X6 
+576X6X^X2 - 18X6X2X1 + 13392X1X2X1 
-144x1X3X5 - 1728X1X3X1^ + 96X2^X3X1 
-216X1X1^X2 - 108X1X5X3^ + 72T^TeT^ 



+774144X1^ + 138240X| + 5308416X2 + 1096704X2 
+1728X2 - 20736Xi^ + 3456X| - 103680Xi^ + 129024X2^ 
+6384X2^ + 1728Xf + 27X3^ + 64X| + 4X| 
+119808X3X1 + 787968X1X1 + 36288X2X32 
+88I28X1X2 + IO8X1X6 + 45888X1X2 
-18432X1^X22 - 9024:T^Tq + 2520T^T^ 
+192X2X2 + 864XiX2 - 432X5X32 - 328X|X6 
-10368X1^X3 - 32X1^X1 + 6592X|X3 - 2976T^T^ 
+432T^T^ + 1728Xf X2 - 432X32X1^ + 1296XfXi 
+864X1^X6 + 432X2X1 - 34560X1^X2 + 32X2^X3 



-IGT^T^ + 224X1X2^ + 4x1X32 



n2 

-1296X5X1X2X3 



If we reexpress these polynomials in the Chevalley basis ri,r2,r3,Te as functions of 
Cartesian coordinates we find 

Pi = -2'' n (3.6) 

A = -2'' ll[sin^{a,x)r (3.7) 

Here we made use of the fact that each orbit Qk,a can be decomposed into a positive 
and a negative semiorbit 



by a hypersurface 



so that 



For example, we define 



(a,0 = (3.9) 



nJ„ = {aGfifc,.;(a,O^0} (3.10) 



^ta = {Aiiei + yU2e2 + yUses + /i4e4 e ^^fc,a; (3.11) 

/ii > or /ii = 0,/i2 > (3.12) 

or /ii = /i2 = 0, /is > (3.13) 

or = /i2 = /i3 = 0, /i4 > 0} (3.14) 

Then the coefficients ^„ of ^ 

e = 5Zae. (3.15) 

n 

must satisfy inequalities, i.e. from Qi and Q2 

> 6 > ^3 > ^4 > (3.16) 

ei>6 + e3 + e4 (3.17) 

^3 implies new inequalities, etc., so that a vector ^ exists for any finite set of orbits. We 
assume ( |3.16| ),( prT^ ) to hold throughout. 

We consider the asymptotic behaviour along an imaginary direction in x-space 

X = -iX^, (3.18) 

In each orbit there exist maximal vectors amax(^, a) so that for fixed ^ 

(amax(/i;,a),0 = max (a,^) (3.19) 

Then (for unique amax(^,a.)) 

Tk,a ~ e^^"-^" (3.20) 

for A 00 (3.21) 



In particular, we have 

amax(l) = ei (3.22) 

amax(2) = ei + 62 (3.23) 
amax(3) = -(361+62 + 63 + 64) (3.24) 
ai^ax(4) = 261 + 62 + 63 (3.25) 

From ( p.l8|) , (|3.19|) we find as leading term in Pi{P2) the unique monomial —T1T^{—T2Tq) 
with asymptotic behaviour 

Pi,2 6^(''^-2'^) (3.26) 

where 

Pl = 56i + 62 + 63 + 64 

P2 = 661 + 462 + 263 (3.27) 

These are special cases of Weyl type vectors 

Pk,a= " (3.28) 

The same leading term ( |3.26| ) results from ( p.6| ), ( p.7| ). In this fashion the constants in 
the factorization formulas ( p.6|) , ( p.7|) can be controlled. 
From (^) follows that 

det g-^ = ^PiP2 = (det J)^ (3.29) 

where J is the Jacobian matrix (|2.22|) . Due to (|3.6|) ,( pr7D det J vanishes in first order 
along the hyperplanes 

(a, x) = 2mT, n E Z,a E Qi^2 (3.30) 

The root space M4 is therefore divided into cells bounded by the walls ( p.30|) where (see 
( p.42|) ) repulsive and impenetrable potentials are positioned. 
Next we derive the r-coefficients from 

^rn — 9mn 'K^ (3.31) 

that ought to be polynomials in the {T„}„g{i_2,3,6}- In fact, we find 

r(i) = (-5Ti - 24, -6T2 - 6T1, -9T3 - I2T2 - 24Ti, -4T1T2 - UTq - 2AT2 + 24Ti) 

(3.32) 

r(2) = (_6Ti, -IOT2 - 48, -I2T3 - 24Ti, -24Tf - 18T6 + 48T3 + 96T2 + 192Ti + 576) 

(3.33) 

From ( p.2| ) and ( |3.32| ),( p.33D we obtain the algebraic differential operator 

« = - E J; + E E '-S'^ (3-34) 

m,n ^ " a=l,2 m ™ 



containing two real coupling constants 7^ as free parameters. This operator is exactly 
solvable in terms of polynomial eigenfunctions 

Dpa(Ti, T2, T3, Te) = Apa(Ti, T2, T3, Tg) (3.35) 

One can show that D possesses a flag of invariant polynomial spaces {V^} 

DVn C Vn (3.36) 

Vn = span{T^'T^^T^'T^<',nieZ>, (3.37) 
ni + 2^2 + 3^3 + 4^6 < iV} (3.38) 

Thus the eigenfunctions {px} can be calculated by linear algebra. 

On the other hand we know from the general scheme that D corresponds by a "gauge 
transformation" ^ to the Schrodinger operator in R4 

H = -A + W (3.39) 

with a Laplacian A and a potential W 

W = \Y.(lalb-\)Rat (3.40) 

a,b 

_1 9 log Pa 5 log 



Rab — ^ g. 



mn 

m,n 



= T^'n^T^^^^n (3.41) 
In Cartesian coordinates we can for P4 express 

Rab = Pab J2 [si^ ^("' ^)]"^ + '^'^b (3.42) 

where 

pn = +1, Cn = -28 

P22 = +2, C22 = -56 (3.43) 

P12 = P21 = 0, C12 = C21 = -36 

can be determined by an asymptotic analysis. 

In this fashion we have derived the Sutherland model for Fa. 
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A Tables 

Table 1: Orbits, generating vectors and cardinals for F4 







#fi 




ei 




24 




ei + 62 


24 


^3 


61 + 62 + 63 


96 




2ei 


, Ci + 69 + e-? + 


24 




2ei 


+ 62 


144 




2ei 


+ 62 + 63 


96 




2ei 


+ 69 + 6q + 6/t 


192 


fig 


2ei 


+ 262 


24 




2ei 


+ 262 + 63 


288 


fig 2 


3ei 




24 


filO 


3ei 


+ 62, 26i + 262 + 63 + 64 


144 


fill 


3ei 


+ 62 + 63 


288 


fil2 


2ei 


+ 262 + 263, 36i +62 + 63 + 64 


96 


^13,1 


3ei 


+ 262 


144 


fil3,2 


2ei 


+ 262 + 263 + 64 


192 


fil4 


3ei 


+ 262 + 63 


192 


fil5 


3ei 


+ 262 + 63 + 64 


576 


fil6 


2ei 


+ 262 + 263 + 264, 461 


24 


^17,1 


3ei 


+ 262 + 263 


288 


fil7,2 


4ei 


+ 62 


144 


fil8,l 


3ei 


+ 362 


24 


fil8,2 


4ei 


+ 62 + 63, 36i + 262 + 263 + 64 


288 


fil9,l 


3ei 


+ 362 + 63 


288 


fil9,2 


4ei 


+ 62 + 63 + 64 


192 


fi20 


4ei 


+ 262, 36i + 362 + 63 + 64 


144 


fi21,l 


4ei 


+ 262 + 63 


576 


fi21,2 


3ei 


+ 262 + 263 + 264 


192 


fi22 


3ei 


+ 362 + 263, 46i + 262 + 63 + 64 


288 


fi23 


3ei 


+ 362 + 263 + 64 


576 


fi24 


4ei 


+ 262 + 263 


96 



Table 2: T„^a expressed as polynomial in Ti, T2, T3, Tg for F4 

T4 = - 2Ts - 6T2 - STi - 24: 
Ts = T2Ti-3r3-6Ti 

TV = TaTi - 472^1 - 8T^ - STq + 22T^ + 36T2 + 8OT1 + 192 

Tg = Ti - 6T^ - 2Tq + 12^3 + 28r2 + 48^1 + 120 
Tg.i = 73T2 - STaTi + 872^1 + 24Ti^ + QTe - 6OT3 - IO8T2 - 228ri - 576 
T9,2 = Tl - ST^Ti - 3T2T1 + STg - 217^3 - 36^2 - QQTi - 192 
Tw = TiT^ - 2T3T2 - ST2T1 - 6Ti - STe - 3OT2 

Tn = TqTi - 2T3T2 + 3T3T1 - 8^271 - 247;^ - 9Te + 63Ts + 108r2 + 240Ti + 576 
T12 = -4^27? - 8Tf - 2TfiTi + 8T3T2 + 12T| + 24T3T1 + 56r2ri + T| + GOT^ 
+40^3 + 120r2 + 288Ti + 288 
^13,1 = T^Ti - 6Tf - 2TeTi - T3T2 + IBT^Ti + I9T2T1 + 2AT^ - 9Te 

+63T3 + IO8T2 + 35IT1 + 576 
^13,2 = - 2T| - TgTi - 6T3T2 - 6T3T1 - 8T2ri - 16r2 - 6Tq 

+2OT3 + 72T2 + 152Ti + 384 
Tu = STsTi^ + 24Ti' + TeTz + GTqTi - 72T3T1 - I36T2T1 - 20r| - 16r3r2 - 3r3^ 

-144T2 + 22T6 - I92T3 - 4OOT2 - 1152Ti - 1536 
= T3T2T1 - 3T3T2 - 4T|Ti - STsT^ - 3T6T2 + TeTi + 4OT3T2 + 21T3ri 

96T2T1 + + 36T| + 24T2 + 9T6 - 6T3 + 84T2 - 240Ti - 576 
Tie = - ATsT^ - 4T2TI + ATqTi + 2T| + 8T3T2 + 6T| - 16T3ri 

-I6T2T1 - 7QT^ + 12T6 - 4OT3 - 72T2 - 4I6T1 - 552 
ri7,i = -2T3T2T1 + 3T3T2 + 4T|Ti + I6T2T2 + 12Ti' + TeTs + 6T6T2 + hT^Ti 

-QTl - 55T3T2 - 72T2^ - 5IT3T1 - 2OOT2T1 - 72Ti^ + 9T6 - I2OT3 

-492T2 - 468Ti - 576 
Ti7,2 = -3T3T2T1 + T2TI - ST^Ti + 3T6T2 - TgTi - 36T| - I9T3T2 

-3T3T1 - 92T2T1 + 24Ti^ + 9Te - 6OT3 - 3OOT2 - 234Ti - 576 
Ti8,i = T| - 2ATf - 15T2T^ - ?>TqT2 - QTqTi + 3T| + 3OT3T2 + 72T3T1 + 54T| 

+ I92T2T1 + 144T2 - 21T6 + I92T3 + 549T2 + 1152Ti + 1536 
^18,2 = nr'l - 2ST2TI - 72Tf - 2nn - 9T6T2 - 26T6T1 + 56T3T2 + 2I6T3T1 + 84T| 

+44OT2T1 + 9r| + 432Ti^ - 84T6 + 576r3 + 13O8T2 + 3456Ti + 4608 
7^19,1 = T3TI - T3T2T1 + 8T2T2 + 4T|Ti - 3T3Tf - 2^^ + STeTs - S^Ti + 6T| 

-8r3r2 - 36T2^ + 1873^1 - 68^2^1 + 487;^ + 18T6 - 69T3 - 408^2 

-480Ti - 1152 

ri9^2 = -8T^ - 4:T2Tf + T^Ti + 9TsT2Ti + 2?>T^Tl + 127^^1 - 2TqTI + 48r2ri^ 

+48rf - TqT^ - ?,TqT2 - ITqTi + 2Tl + 24T3T2 + 367^ + 7273^1 + 244T2T1 
+384Ti2 + I8T3 + I92T2 + 512Ti 
T20 = -17167^2 + 48r6-2r3r|-8ri7f-2ri2r6 - 29287^1 - 636^3 

+76Ti^T2 + T^Ti - 22Kri - I52T2T3 + 22T{rQ - 2UT]Tz - 306Ti^ 
-68071^2 + 127^2^6 - 21Tl + 2ATIT^ + 1207f + 473^6 - 67f - 3600 - 



= -4608 - 1008T2 + lllTg - 576T3 - 3456Ti + UTiT^^ + AT^T^^ + STsTl - IGTsTa 
-2I6T1T3 - I52T1T2 + 88TfT2 + 56TiT| - 2T'ITq - 8TfT2 + ISTaTg + TsTI 
-AT^Tl - 2TqTiT2 + 24T1T2T3 - 432Ti2 - 9X3^ + 36r| + 72ri^ + 12r| 

= 4320 + IIO4T2 - 152T6 + 912T3 + 3648ri - 64Tir6 
-I6T3T6 - 8X3X2^ + 96T2T3 + 384TiT3 
+I28T1T2 - 208T^T2 - 80TiT| + ST^Tq 

+i6T^T2 - 24T2T6 - 2T2TI + at^t:^ 

-AST^Ts + 4TqTiT2 - 48T1T2T3 + 312X^2 
+48T| + 12T| - 192Tf - 8T| + 12T^ + T| 



